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ABSTRACT 

We demonstrate that a magneto-convection simulation incorporating essential 
physical processes governing solar surface convection exhibits turbulent small- 
scale dynamo action. By presenting a derivation of the energy balance equation 
and transfer functions for compressible magnetohydrodynamics (MHD), we quan- 
tify the source of magnetic energy on a scale-by-scale basis. We rule out the two 
alternative mechanisms for the generation of small-scale magnetic field in the 
simulations: the tangling of magnetic field lines associated with the turbulent 
cascade and Alfvenization of small-scale velocity fluctuations ( "turbulent induc- 
tion"). Instead, we find the dominant source of small-scale magnetic energy 
is stretching by inertial-range fluid motions of small-scale magnetic field lines 
against the magnetic tension force to produce (against Ohmic dissipation) more 
small-scale magnetic field. The scales involved become smaller with increasing 
Reynolds number, which identifies the dynamo as a small-scale turbulent dy- 
namo. 

Subject headings: Sun: magnetic fields - MHD - turbulence - Sun: photosphere 



Introduction 



Generically, three-dimensional (3D) magnetohydrodynamic turbulence excites dynamo 
action when the magnetic Reynolds number exceeds a critical threshold (such that ampli- 
fication by stretching dominates over Ohmic dissipation). That turbule nce could amplify 
magnetic energy by the random stretching of field lines wa s proposed by iBatchelorl (119501 ) 
and first demonstrated in direct numerical simulations by Meneguzzi et al. (jl981 ). Many 
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turbulent flows with high enough Reynolds numbers support dynamo action. Solar sur- 



face convection is 



clear (IQssendriiver 



ikely such a flow, but the relation to the global solar dynamo is not 



2003J). There is evidence for local dynamo action near the solar surface 



( jPetrovay Sz Szakalyl Il993l ). and the question of whether surface convection can support a 
turbulent dynamo is the focus of the present work. 



The simulations of iMeneguzzi et al.l ( 119811 1 demonstrated turbulent dynamos with dra- 
matically different characteristics. In one case, magnetic energy grows at scales smaller than 
the forcing scale of flu id motions. This de fines small-scale dynamo (SSD) or fluctuation dy- 
namo action (see, e.g.. llskakov et al.ll2007l ). In the other case, magnetic energy grows at scales 
larger than the forcing scale: large-scale dynamo. Such large-scale dynamos are often studied 
in the framework of mean-field theory (se e, e.g.. lKrause fc Raedlerlll980l ; ICovas et al.lll997l ; 
Field et al.lll999l ; iField fc Blackmanl 120021 ) which suggests that the production of large-scale 
magnetic energy is related to the lack of reflectional symmetry in the small-scale velocity 
fluctuations (e.g., from helical motions). A more intuitive picture is that kinetic helicity cre- 
ates magnet ic helicity through A lfvenization, which then goes through an inverse cascade to 
large scales ( Pouquet et al.|[l976l ) - the magnetic tension force expands coils in the magnetic 
field lines. Symmetry breaking (e.g., helicity) is not required for small-scales dynamos and it 
is widely believed that sufficiently chaotic 3D flows will be SSDs. Near the solar surface, the 
convective (granulation) time scale is much shorter than the rotation period, so that a flow 
with negligible net helicity results and large-scale dynamo action is not expected. However, 
surface convection is expected to be sufficiently chaotic to be a small-scale dynamo. 

The origin and properties of the quiet-Sun intra-network magnetic field provide obser- 
vational evidence of and argu ments for local small - scale dynamo action on the Sun. Firstly, 
from a flux-transport model, iPetrovay &: Szakalyl (119931 ) conclude that the decay of active 
regions (or of flux tubes not quite rising to the surface) is insufficient to account for the 
observed intra-network magnetic fields. A source term is required in their model to match 
the observations and they conclude that the physical interpretation of this term is small-scale 
dynamo action in the convection zone. Secondly, small-scale dynamo action is also consis- 
tent with high resolution magnetograms of the intra-network. These show mixed polarity 
fields on small scales (variously called the m agnetic carpet or the salt-and-pepper pattern; 
Title fc Schrijverlll998l ; iHagenaar et al.l 120031 ). The fractal n ature of these opposite polar i- 



ties extends down to the resolution limit of the observations ( jPietarila Graham et al. 



2009|). 



Furthermore, the amount of observed small-s cale flux is not dependent on the solar cycle 



nor does it show any latit udinal dependence ( IHagenaar et al.l 120031 ; I Sanchez Almeida] 12003 



Trujillo Bueno et al.ll2004j ). These results all suggest that the source of quiet-Sun magnetic 
field is independent of the global dynamo. Lastly, turbulent convection has been shown 
to drive dynamo action in numerical simulations of Boussinesq convection without rota- 
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tion (jCattanedll999l ; ICattaneo et al.ll2003l ). Observations and simulations together provide 
evidence, then, of a small-scale dynamo driven by turbulence at the solar surface^ 

An alternative interpretation of the observations is that the small-scale field results 
from turbulence acting on the large-scale magnetic field from the global dynamo. Such 
induced small-scale fields could result from two different processes. One process is the tur- 
bulent tangling, sometimes called "shredding," of field lines which moves magnetic energy 
to smaller scales as part of the turbulent energy cascade. However, the turbulent cascade 
is associated with power-law energy spectra and, therefore, the amount of small-scale intra- 
network flux should change as the strength of any large-scale backg r ound field changes over 



the solar cycle. This contr adicts observations ( IHagenaar et al.ll2003 



Sanchez Almeida 2003 



Trujillo Bueno et al.l 120041 ) . A similar argument, agai nst the decay of active r egion s as the 
source of intra-network flux, has been put forward by lSanchez Almeida et al.l ( 120031 ). Addi- 
tionally, the total unsigned magnetic flux (and energy) in active regions even during solar 
maximum is less than the unsigned magnetic flux (and energy) contained in the quiet-Sun 
( iSanchez Almeidal 12004 ; iTrujillo Bueno et al.l 120041 ). making decay from active regions as a 
source of the small-scale field very unlikely. 

The second proc ess for producing small-s cale magnetic field from large-scale field, called 
turbulent induction (jSchekochihin et al.ll2007l ). involves the stretching of a uniform (or large 
scale) background magnetic field by turbulent fluid motions which excites Alfvenic magnetic 
fluctuations on the same scale as the turbulent fluid eddies. Alfvenic turbulent induction 
will be present whenever there is a significant background field. In the presence of both a 
flow capable of sustaining small-scale dynamo action and a large-scale magnetic field@ the 
small-scale dynamo and Alfvenic induction may compete as the source of small-scale field. 

Given the prevalence of turbulent dynamos in homogeneous, isotropic turbulence, such 
dynamo action is expected in the Sun unless additional physics can be identified which would 
suppress it. Two points are often raised to argue that a small-scale dynamo cannot operate 
in the Sun. Firstly, turbulent eddies smaller than the characteristic scale of the magnetic 
field act like a turbulent magnetic diffusivity and could inhibit dynamo action. This is 
an important concern for the Sun (and other stars) because the kinetic Reynolds number, 
Re = VqIq/v (i>o and Iq being typical velocity and length scales, and v the kinematic viscosity), 
is much larger than the magnetic Reynolds number, Re u = vqIo/v (v being the magnetic 



1 Actually, simu lations suggest turbulent dynamo action occurs in the bulk of the convection zone as well 
(jBrun et alll2004h - 



2 For very strong background fields (several times the eq uipartition field strength, B ea m y/4npoV r 
350 G), the large-scale magnetic field quenches the dynamo (jHaugen fc Brandenburg||2004l ). 
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diffusivity) . Their ratio, the magnetic Prandtl number, Pm = Rcm/RGi is approximately 
10 -5 near the solar surface. The critical magnetic Reynolds number for turbulent dynamo 
action, Re^, sharply increases with decreasing P M and it has been suggested that Re M 
goes to infinity as Pm goes to zero fjRogachevskii Sz Kleeorinl 119971 ; iBoldyrev fc Cattaneo 
20041 ; ISchekochihin et al.ll2005l ). However, recent numerical simulations a ttest that Re^ ap - 
proaches a plateau for Pm -C 1 both with ph ysical (Laplacian) viscosity (jPonty et al.l 120051 ) 
and with hyperviscosity (jlskakov et al.ll2007l ). Baring the identification of a new length scale 
in the problem, these results indicate that, for magnetohydrodynamics (MHD), small-scale 
dynamo action remains possible in the asymptotic limit as Pm goes to zero. Additionally, 
a laboratory dynamo with liq uid sodium (Pm y 10~ 5 ) resulting from unconstrained turbu- 
lence has been demonstrated (IMonchaux et al.l 120071 ). This establishes that a turbulent (at 
l east, a large-scale) dy namo is possible at values of Pm corresponding to the solar plasma 
f lMonchaux et al.ll2009| ). 



The second suppression argument is that, unlike the lCattaneol ( 119991 ) Boussinesq simula- 
tion, the Sun is strongly strat ified and magneti c flux is swept into the downflows and subject 
to long reci rculation times (|Stein et al.ll2003l ). In their stratified simulations with open 
boundaries, IStein et al.l (120031 ) found no evidence of dynamo action for surface convection. 
Their magnetic Reynolds number, however, was below the critical value. IVogler fc Schiissler 
( 120071 ) have demonstrated that surface convection with little field recirculation and strong 
density stratification (as well as other physical effects present in the Sun) can support local 
dynamo action when Rcm ^ 2000. We will determine if this local dynamo action is a small- 
scale turbulent dynamo. Currently, no other likely suppression mechanisms for small-scale 
dynamo action in the intra-network photosphere are known. 



The magnetic energy spectrum of the lVogler &: Schusslerl (120071 ) dynamo peaks at scales 
smaller than the energy-containing scale of the fluid motions, which is suggestive of a small- 
scale dynamo. However, except in the most idealized of staulationsg non-zero time-averaged 
mean flows exist for times (~ 10 minutes for photospheric convection) long compared to 
inertial-range eddy-turnover times (e.g., ~ 6 seconds at a scale of 1km). Therefore, even 
in the absence of net helicity, this mean flow can act as a low ReM dynamo that pro- 
duces magnetic field near the mean-flow scale (jPonty et al.l 120071 ) , ~ 1 Mm for convective 
granulation. This would not be a small-scale dynamo; it would still have small-scale mag- 
netic field produced from either the turbulent cascade or from Alfvenic turbulent induction. 
In order to fully understand what is occurring, it is important to disentangle the possi- 
ble sources of small-scale magnetic energy and properly identify the dynamo mechanism 



3 TIie case of delta-correlated in time, isotropic, and homogeneous forcing being the sole exception. 
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( ISchekochihin et all 120071 ). 

We have employed transfer analysis to measure, scale by scale, the relative strengths 
of the sources of magnetic energy: turbulent cascade, Alfveni c turbulent induction, an d 
dynamo action (by stretching of field lines) in the simulations of IVogler fc Schiisslerl ( 120071 ). 
The scales involved in the dynamo and their dependence on Reynolds number as well as 
growth rates and energy spectra allow us to determine that the dynamo is a turbulent small- 
scale dynamo. 



2. Data and methods 



We use MURaM (jVoglerll2003l ; IVogler et al.ll2005l ) to perform simulations for a rectangular 
domain of horizontal extent 4.86x4.86 Mm 2 and a depth of 1.4 Mm (800 km below and 600 km 
above the simulated solar surface). The sides boundaries are periodic in both horizontal 
directions. The open lower boundary assumes upflows to be vertical, v x = v y = d z v z = 0; for 
downflows vertical gradients are set to zero, d z v x = d z v y = d z v z = 0. The upper boundary is 
closed. The magnetic field is vertical at upper and lower boundaries: B x = B y = d z B z = 0. 
This ensures that there is no Poynting flux into or out of the box. The magnetic diffusivity is 
i ncreased in the lower 150 km of the box in order to well resolve the diffusive boundary layer 
( IVogler fc Schiisslerl 120071 ) . These boundary conditions allow the convective downflows to 
leave the box and, thus, simulate an artificially isolated surface layer. These same downflows 
drag magnetic field to the lower region of enhanced diffusivity where it can be eliminated 
(this simulates the fact that the field should not be available for further stretching due to 
long recirculation times; see lStein et al.ll2003l ). The upper 600 km of the box is convectively 
stable and the lower 800 km of the box is convectively unstable. The convection is driven by 
radiative cooling (calculated using grey radiative transfer and the Rosseland mean opacities) 
at the surface where the optical depth is unity. The effects of partial ionization are included 
in the equation of state. These simulations are then so-called "realistic" simulations of 
a portion of the solar surface layer including compressibility and stratification (4 orders- 
of-magnitude variation in the density), radiative energy transport, partial ionization, and 
little recirculation. "Realistic" is used here to distinguish them from a class of "idealistic" 
simulations of incompressible, isotropic, homogeneous, and triply-periodic MHD turbulence. 
Of course, neither type of simulation is able to achieve the Reynolds and Prandtl numbers 
of the Sun, but aside from this, MURaM simulations include more of the physics relevant to 
the near-surface layers of the Sun. Our aim is to determine if these effects inhibit turbulent 
small-scale dynamo action or overshadow it with some other mechanism. 



Dynamo runs with increasing resolution and Reynolds numbers have been carried out. 
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Table [T] summarizes the parameters of the runs. Figure [T] displays the time evolution of 
m agnetic energies. R u n 2 h as previously been reported (as "Run C") for dynamo action 



in 



Vogler fc Schusslerl (120071 ). All runs start from an initial hydrodynamic case plus weak 



magnetic seed field. There is an initial growth of Em for about 3 minutes due to flux expulsion 
and convective intensification followed by the linear (kinematic) phase during which magnetic 
energy grows exponentially with time and the magnetic field is too weak to affect fluid 
motions. The lower resolution simulations are continued until the nonlinear phase when the 
back reaction of the Lorentz force on fluid motions becomes important and the magnetic 
energy begins to saturate. 
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Fig. 1. — Mean magnetic energy density, (e^), versus time, t, for dynamo simulations. 
Shown are Run 1 (red, dotted), Run 2 (black, solid), Run 3 (blue, dashed), and Run 3-P 
(cyan, dash-dotted). 



For all MURaM simulations, constant magnetic diffusivity is employed (outside the re- 
gion of enhanced diffusivity near the bottom boundary), bu t artificial shock re solving; and 
hyperviscosity are applied to the momentum equations (see IVogler et al.ll2005l for details). 
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This approach allows us to reach high effective kinetic Reynolds number, Re c f[, without a 
prohibitively small time step. No value of viscosity is defined, but Re e s can be estimated 
from the energy spectra (Figure The estimate of the effective Reynolds number is derived 
from the representative length scales: the velocity Taylor microscale, \y, 



XI = (, 2 >/<u; 2 > 



roo roo 

' E v (k ± )dk ± / / klE v (k ± )dkj 
o Jo 



(1) 



where Ey(k±) is the horizontal velocity spectrum@ and the integral scale for the turbulent 
motions, L , 

poo roo 

L = k- l E v (k ± )dkJ E v (k ± )dk ± . (2) 



The e ffective Reynolds numbers of the simulations is given by (jBatchelorlll953l ; IWeygand et al 
2007h . 

' T, o 



Re P w oc 



(3) 



where the constant of proportionality is unknown. We can also determine the magnetic 
energy Taylor scale, \ M , 



E M (k ± )dk ± / / k 2 ± E M {k ± )dk A 



\ 2 

10 Jo 
This allows us to estimate the effective magnetic Prandtl number, 



PM,eS OC — . 



(4) 



(5) 



We determine the constant of proportionality from an incompress ible, isotropic, and h omoge- 
neous dynamo simulation. We employed a pseudospectral code ( iGomez et al.ll2005al Jbl) and 
forced the velocity field with a superposition of harmonic modes with random phases, the 
resolution was 256 3 grid points, and Pm = 4. From this calibration of Eq. we estimated 
that the magnetic Prandtl numbers in our MURaM runs are between 1 and 2 (see Tabled]). 

We find power- law scalings, kj_, in the inertial range with (3 = — 1.2 ± 0.1 for the 
kinetic energy spectrum and /3 = —1.54 ± 0.07 for the velocity spectrum. Using second- 
order longitudinal structure functions for the velocity (not shown), we determined that 
at scales larger than w 50 km (approximately half a pressure scale height of the vertical 
stratification at the surface) the flow is anisotropic. The strong vertical downflows in the 



We employ horizontal spectra which are obtained by performing a two-dimensional Fourier transform 
of any field, /, at each horizontal layer resulting in f{k x ,k y ,z,t). This quantity is then projected onto a 
one-dimensional wavenumber k\ = k 2 x + ky. The horizontal spatial frequency, v±, is k±/2ir. 
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Table 1. Properties of MURaM dynamo runs: grid resolution, effective magnetic Prandtl 
number, Pm, c s, magnetic diffusivity, rj, magnetic Reynolds number, ReM based on 
v rms ~ 3.3km/s and a length scale of 1 Mm, and kinematic dynamo regime e— folding time, 



Simulation 


Grid Resolution (km) 


PM,e« 


7] (cm 2 s x ) 


Re M 


7- 1 (s) 


Run 1 


9 x 9 x 10 


~ 2.0 


1.6 ■ 10 10 


« 2100 


1600 


Run 2 


7.5 x 7.5 x 10 


~ 1.3 


1.25 ■ 10 10 


« 2600 


660 


Run 3-P 


5x5x7 


~ 0.8 


1.25 • 10 10 


« 2600 


1200 


Run 3 


5x5x7 


~ 1.1 


6.25 • 10 9 


« 5300 


230 


Run 4 


4x4x4 


~ 1.1 


4 • 10 9 


« 8300 


150 



km km 
10000 1000 100 10 10000 1000 100 10 




10~ 9 10~ 8 10~ 7 10~ 6 10~ 9 10~ 8 10" 7 10" 6 

Uy (cm -1 ) Vi (cm" 1 ) 



Fig. 2. — Left: Kinetic energy spectrum, Ek, versus horizontal spatial frequency, u±, for 
Run 2 (black, solid) averaged over t € [2.2, 2.4] hr, for Run 1 (red, dotted) averaged over 
t G [0.2, 1.4] hr, and for Run 3 (blue, dashed) averaged over t G [0.1, 0.5] hr. The vertical 
lines represent the integral scale (~ 10 8 cm) and Taylor scales near the beginning and end 
of the inertial range where we identify a power-law scaling: /cj 1 ' 2 (dot-dashed line). Right: 
Velocity spectrum, Ey, versus u± for the same runs. The identified power law extends to 
smaller scales and has a slope of —1.54 ± 0.07. 
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convectively unstable lower 800 km of the box are the source of this anisotropy. Because 
of the anisotropy, our spectra cannot be compared to th e = —5/3 from the Kolmogorov 
spectrum for homogeneous, isotropic turbulence (see, e.g. jFrischlll995l ). 



3. Comparison with "idealistic" small-scale dynamos 
3.1. Incompressible, isotropic kinematic/linear phase 



In the kinematic regime, where the Lorentz force is insignificant, the magnetic energy is 
observed to grow exponentially, Em ~ exp7t (7 is the growth rate). For turbulent dynamos 
(with Pm Cl), the growth rate scales as 

1/2 



7 



Re 



M 



(6) 



thou gh no simulation has yet exceeded Re M enough to observe the predicted scaling (see, 
e.g. 



Schekochihin et al.l 120071 ) . For the large magnetic Prandtl case, Pm > 1, we expect 



Re 1 ^ 2 . Again, this scaling has not yet been found in any simulations ( Schekochihin et al. 



20041 ). Other predictions about the kinematic (line ar) phase of the small-sca le dynamo are 
due to the exactly solvable Kazantsev model (see I Schekochihin et al.l 120041 and references 
therein). This model predict s that the amplitude o f each (Fourier) mode grows exponentially 
at the same rate (see also, iMininni et al.l l2005bl for a growth rate analysis) and that the 
magnetic energy spectrum at large scales follows a k 3 ^ 2 scaling. 



3.2. MURaM kinematic/linear phase 

For the MURaM dynamo, the dependence of growth rate on Reynolds number is shown in 
Figure [3j The growth rate is an increasing function of Reynolds number. This indicates that 
the dynamo is a smal l-scale (inertial range) pro cess. We find that we have not yet sufficiently 



exceeded Re M (from lVogler fc Schusslerl 120071 and Run 1, 1300 < Re M < 2100) to observe 



the predicted power-law scaling. Instead, at the Reynolds numbers we can achieve, the 
growth rate depends almost linearly on ReM (dashed line). Higher resolution simulations 
will be required to determine if there is, indeed, an asymptotic power-law. 

In Figure HI we have plotted magnetic energy spectra, E M (v±), at different times from 
the beginning of the linear regime (including the end of the flux expulsion phase) for Run 3. 
Power laws can be fit to the magnetic spectrum for scales larger than Xk, where magnetic 
energy is growing at the same rate for all scales. A power law is in contrast to the case of 
a laminar dynamo and indicates that a turbulent, self-similar process underlies the dynamo 
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Fig. 3. — Growth rate of magnetic energy, 7, versus magnetic Reynolds number, Rcm for 
Run 1, Run 2, Run 3, and Run 4. Also shown is a linear fit (dashed) and the theoretical 
Re]^ 2 scaling (dotted) for reference. 
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mechanism. We find EM{k) oc kj_ with (3 = 0.83 ± 0.04. The results from Run 2 are 
(3 = 0.63 ± 0.06 and for Run 1, (3 = 0.53 ± 0.05. None of these slopes agree with k 3 ^ 2 
from the isotropic Kazantsev case, but are not expected to because of the anisotropy of 
fluid motions in the downflows. We find, however, that the degree of anisotropy decreases 
with increasing Re e g for the runs and the exponent becomes closer to 3/2. This decrease 
in anisotropy is due to the generation of stronger horizontal fluctuations by the turbulence 
(strong vertical fluctuations are induced by the convective driving). In summary, we find 
that the general character of the MURaM dynamo is similar to the incompressible, isotropic, 
small-scale dynamo but the spectral index differs. 
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Fig. 4. — Magnetic energy spectra, Em{v±) (solid lines), for Run 3 for t = 11, 16, 23, 28, 
and 31min (kinematic phase). The vertical dashed line is the kinetic Taylor scale. Above 
this scale, the magnetic energy spectra follow a power law (indicating a turbulent dynamo 
process) with all scales growing at the same temporal rate. The time-averaged kinetic energy 
spectrum (dashed line) is compensated by k]^ 2 so the inertial range is flat. 



4. Transfer analysis (kinematic phase) 



It is clear that the magnetic energy is peaking at scales nearly one order of magni- 
tude smaller than the energy- containing scale of the fluid motions, Lo ~ 1 Mm (see Figure 
H|). This indicates small-scale dynamo action, with one caveat. As discussed in Section 
[U "mean flows" might give rise to a Mm-scale dynamo. We must, therefore, disentangle 
three possible sources of small-scale magnetic energy: the tangling of field lines caused by 
the turbulent energy cascade, Alfvenic response of larger-scale field to smaller scale velocity 
fluctuations (Alfvenic tur bulent induction), and dynamo stretching of magnetic field lines 
(ISchekochihin et al.l 120071 ) . This can be accomplished by spectral transfer analysis. 



Spectral transfer analysis was introduced by iKraichnanl (119671) and is widely used to 
understand incompressible turbulent p rocesses in both two (e.g 



Maltrud fc Vallis 



1971 



Evinkil2006[) and three (e.g ., IZhoulll993l ) dim ensions for both Navier-Stokes (e.g.. 



Mininni et all 120081 ) and MHD (e.g., iDebliquv et al.l 120051 : IVerma et al. 



1993 



Kraichnan 



20051 ). In 



Appendix IA.ll we extend it to compressible MHD. Transfer analysis allows us not only to 
quantify the sources of magnetic energy but also the scales at which they operate and the 
scales at which they generate magnetic energy. In general, the transfer function TxYF{k) 
measures the net rate of energy transfer from energy reservoir X to reservoir Y mediated by 
force F. For T X YF(k) > (< 0) net energy is received (lost) by Y at the scale given by 2ir/k. 
Transfers between the kinetic and magnetic energy reservoirs occur via the Lorentz force, 
which can be separated into the effects of a magnetic pressure and a magnetic tension force. 
For example, T B Kr{k) measures the net work done by the magnetic tension force on the 
fluid motions at wavenumber k by all scales of the magnetic field. When it is negative, net 
kinetic energy is lost by fluid motions at wavenumber k working against the magnetic tension 
force. Likewise, Tx B r{k) measures the net work done on the magnetic field at wavenumber 
k by all scales of fluid motions. These two transfer functions measure energy exchanged by 
the kinetic and magnetic energy reservoirs and, therefore, the rate of energy gained by the 
magnetic field is equal to the negative of the rate of energy lost by fluid motions, 



^2 T KBr{k) ^ T BKT (k) 



(7) 



Similarly, —T BKP {k) measures work against magnetic pressure gradients and T KBP {k) the 
net magnetic energy generated, with 



}T K Bp{k) = - T BKP {k) 



This second transfer is peculiar to compressible MHD as J^ fc TxBp{k) = for the incompress- 
ible case where T^bp measures only the effects of the turbulent cascade moving magnetic 
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energy to smaller scales. In our case Tkbp measures both the magnetic portion of the 
turbulent cascade and magnetic energy generated by compressive motions. 
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Fig. 5. — Net energy transfers in the kinematic phase. Work against magnetic tension, 
Tbkt < (stretching), accounts for 95% of magnetic energy generated while work against 
magnetic pressure, Tbkp < (compression), accounts for 5%. The latter is involved in 
the process of breaking down larger-scale field, Tkbp < (25%), into smaller-scale field, 
Trbp > (30%), as part of the turbulent cascade. The dominant producer of magnetic 
energy is the stretching of magnetic field lines against the magnetic tension force, Tkbt > 
(associated with turbulent dynamo action). 

In Figure [51 we illustrate the net energy transfers between kinetic and magnetic energies, 
and in Figure [6] we plot their scale dependence. We find that the dominant source (95%) 
of magnetic energy is the stretching of field lines against the magnetic tension force. A 
lesser source of magnetic energy (5%) is work against the magnetic pressure force. This 
latter energy is combined with larger-scale (> 30 km) magnetic energy lost (25%) to the 
"magnetic cascade" to produce (30%) smaller-scale (< 30 km) magnetic energy. Using a 
similar expression to Equation (T5]), we determine that the predominant scale at which fluid 
motions are doing work against the magnetic tension force to stretch magnetic field lines is 
140 km. As seen in Figure E] for Run 3, —T BK t peaks at the corresponding spatial frequency, 
7-10 -8 cm" 1 . Acceleration of large-scale fluid motions by the magnetic tension force, Tbkt > 
green dash-dotted line, shown in Figure [6] is mostly likely due to the somewhat artificial 
separation of the Lorentz force into magnetic tension and an isotropic magnetic pressure. 
As there is no Lorentz force along magnetic field lines, there can be no transfer along field 
lines and a portion of negative Tbkp is offset by positive Tbkt- To see at which scales the 
magnetic field is gaining or losing energy, we examine Tkbt- We find that net magnetic 
energy is gained at all scales of the simulation. The predominant scale for magnetic energy 
production by stretching is identified as 65 km (a spatial frequency of 1.5 • 10~ 7 cm _1 ). Work 
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against magnetic pressure, Tbkp < shown as a cyan dotted line, is mainly by granulation- 
scale fluid motions. The net result is to remove (T^bp < 0, red dotted line) larger-scale 
magnetic energy, which was produced by stretching motions, and break it down into smaller- 
scale magnetic structures {T^bp > 0, black solid line). 
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Fig. 6. — Net energy transfer rates for Run 3 versus horizontal spatial frequency, v±, av- 
eraged over t G [ll,28]min (linear dynamo regime). Left: Work against magnetic ten- 
sion, TsKT{k) < 0, is shown as pink dashed lines; fluid acceleration by magnetic tension, 
TsKT^k) > 0, as green dot-dashed; and work against magnetic pressure, T B Kp{k) < 0, as 
cyan dotted. Right: Dynamo stretching, Tkbt > 0, is shown as blue dot-dashed lines; 
generation by compression (and received from other scales), T^bp > 0, as black solid; and 
magnetic energy removed by compression, T^bp < as red dotted. The error level of the 
transfer analysis is 2.5 • 10 2 erg cm~ 2 s _1 (see Appendix IA.2I) . 



With the sources of magnetic energy quantified on a scale-by- scale basis, we can now 
identify the source of the small-scale magnetic energy seen in Figure HI This energy is seen to 
peak at scales between 60 and 120 km. As we can see in Figure [6J the "cascade" dominates 
the production of magnetic energy only at scales below 20 km. Just as for incompressible 



turbulent dynamos (IMininni et al.ll2005al ). there exists a range of scales where the amplifica- 
tion of the magnetic field is dominated by injection of energy from turbulent stretching while 
the magnetic cascade dominates only at smaller scales. The source of the ~ lOOkm-scale 
magnetic energy is a dynamo process and not the turbulent tangling of larger-scale field 
lines. As the magnetic energy is predominantly produced at a scale of 65 km due to the 
stretching by fluid motions which operate predominately at a scale of 140 km (the peak of 
—Tbkt in Figure [6]), this suggests that the field lines stretched by these motions have a scale 
between 45 and 120 km. This is due to the fact that all transfers occur between three fields 
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whose wave-vectors must form a triad (see Figure 

ED 13 This allows us to rule out Alfvenic 
turbulent induction, i.e., small-scale velocity fluctuations interacting with a large-scale field 
to produce small-scale magnetic energy, as an important source of the small-scale magnetic 
energy in the simulations. Furthermore, as the three scales involved in the dynamo all lie in 
the inertial range, we can identify the MURaM dynamo as a turbulent small-scale dynamo. 




Fig. 7. — Illustration of triadic transfers. Left: Case for the linear phase of the MURaM 
dynamo Run 3. Fluid motions predominately at a scale of ~ 140 km create magnetic energy 
predominately at a scale of ~ 65 km. As the three wave- vectors must form a triad, the scale 
of the magnetic field being stretched must have a scale of 80 ± 40 km. Right: The Alfvenic 
case: small-scale velocity fluctuations interact with a large-scale field to produce small-scale 
magnetic energy. 

For a small-scale dynamo, the smallest eddies provide the fastest stretching and this 
should be reflected in the transfer analysis. We would expect both the dominant scale of 
stretching motions and of magnetic field production to decrease as the Reynolds numbers 
are increased. That is, we expect the dynamo action to move to smaller and smaller scales. 
This is, in fact, what we find in the analysis. The dominant scale of stretching motions is 
200 km for Run 1, 180 km for Run 2, and 140 km for Run 3 while the dominant scales of 
magnetic field production are 110 km, 100 km, and 65 km. This provides further evidence 
that the MURaM dynamo is a turbulent small-scale dynamo. 



5 See Equation (|A9[) . describing the production of magnetic energy at wavenumber k. The expression 
involves the transform of the product of the magnetic field (at some scale, p) and the velocity field (at 
another scale, q). By the convolution theorem, we see that p must be given by k — q. 
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5. Nonlinear /saturated phase 

The analysis of the nonlinear phase of the dynamo is quite similar to that already 
presented for the linear phase. For the sake of brevity, we describe only the differences 
between the two phases without figures. The magnetic energy spectra peak at scales slightly 
larger than the kinetic Taylor scales, \k- For scales larger than A#, a power-law kj_ with 
f3 = 0.4 ± 0.1 can be fit. This is not as steep as in the linear regime and does not appear 
to be sensitive to the Reynolds number. It might be sensitive, however, to the large-scale 
flow geometry which will be investigated in a future work. There is a net transfer of energy, 
at 2% for Run 1 (4% for Run 2 and 8% for Run 3) of the total magnetic energy generation 
rate, to fluid motions via the magnetic tension force for scales smaller than 50 km. These are 
short-wavelength Alfven waves. Magnetosonic waves are also produced, at 0.5% for Run 1 
(1% for Run 2 and 2% for Run 3) of the magnetic energy generation rate, for scales smaller 
than 100 km. The saturation of the large-scale magnetic energy is reflected by a near balance 
of generation by stretching and losses to the compressive "cascade," T KBT {k) + T KBP {k) ~ 
for wavenumbers corresponding to scales above a few hundred km. We also note that the 
representative scale at which fluid motions are doing work against the magnetic tension 
force has increased, for Run 2, from the linear phase (where it is 180 km) up to 250 km. 
This result is similar to the result for an incompressible small-scale dynamo where in the 
saturated phase (as opposed to the linear p hase), forcing-scale e ddies dominate the energy 



injection as compared to turbulent eddies (lAlexakis et al.ll2005l ). The representative scale 



for the production of magnetic energy in Run 2 is 130 km in the saturated state. From 
triadic considerations, we expect that predominately magnetic field with scales between 85 
and 270 km is stretched. All results are again indicative of a small-scale turbulent dynamo 
as all scales in the triad are much smaller than the ~ 1 Mm energy-containing scale of the 
granulation convection. 



6. P M < 1 

While our dynamo calculations include much of the physics believed to be present in the 
solar photosphere, the limited computational resources available today do not accommodate 
realistic values of the magnetic and kinetic Reynolds numbers. Moreover, the solar value of 
their ratio, the magnetic Prandtl number ~ 10~ 5 , will remain unachievable for strongly strat- 
ified, radiative small-scale dynamo simulations for many years to come. We can, however, 
take the first steps in this direction by employing the grid size of Run 3 with the magnetic 
diffusivity of Run 2. This run, Run 3-P (initialized from 15 minutes into Run 3 to reduce 
the computational expense) has Pm ~ 0.8 and an exponential growth rate of 7 ~ 1/1200 s^ 1 
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(see Table [T] and Figure [T]). The lower growth rate compared to Run 2 is likely due to a 
decre ased ratio of ReujR^M (R e M increases with decreasing P M , e.g., ISchekochihin et al. 



20051 ) . The results of the transfer analysis are very similar to those shown in Figure El We 



identify the dominant scale of stretching motions as 150 km (close to that of Run 3) while 
the dominant scales of magnetic field production is 85 km (intermediary to Run 2 and Run 
3). This indicates that the magnetic field being stretched has a scale of 130 ±70 km. We con- 
clude that the small-scale turbulent dynamo action is essentially unchanged in comparison 
to the runs with P M > I. 



7. Summary and Conclusions 

We have determined that the MURaM surface dynamo is a turbulent small-scale dynamo. 
The dynamo growth rate, 7, increases with Reynolds number, consistent with the picture of 
dynamo action moving to smaller scales. At scales larger than the peak of magnetic energy, 
the magnetic spectrum is a power-law indicating a self-similar (e.g., turbulent) dynamo 
mechanism. The magnetic energy spectrum peaks at scales an order of magnitude smaller 
than the energy-containing scale of fluid motions (the granulation scale). This peak moves 
to smaller scales with increasing Rcm- By deriving spectral transfer analysis for anisotropic, 
compressible MHD, we are able to identify the source of this small-scale magnetic energy. 
The analysis ruled out two other possible mechanisms for the generation of the small-scale 
magnetic field: the tangling of larger-scale magnetic field lines associated with the turbulent 
cascade and Alfvenization of small-scale velocity fluctuations (or "turbulent induction"). We 
demonstrated that the source of the magnetic field was due, rather, to stretching motions in 
the inertial range of the turbulence. The inertial motions stretch predominantly small-scale 
magnetic field to produce more small-scale magnetic field. All three scales are significantly 
smaller than the granulation scale and move to even smaller scales with increasing Reynolds 
numbers. This positively identifies the dynamo process as a small-scale turbulent dynamo. 

We also identified the key differences between the MURaM small-scale dynamo and those 
studied in the isotropic, incompressible MHD case. The presence of compressive motions 
opens up a new mechanism for the transfer of kinetic to magnetic energy. Namely, compres- 
sion becomes involved in the "cascade" of magnetic energy to smaller scales, resulting in a 
net increase (summed over all scales) of magnetic energy. This mechanism accounts for only 
5% of the rate of magnetic energy generated (and slightly decreasing with increasing Rcm)- 
It is accomplished by Mm-scale fluid motions compressing larger scale magnetic structures 
(which are generated by stretching motions) into scales much smaller (< 20 km) than the 
peak in the magnetic energy spectrum. This magnetic energy is then, presumably, dissipated 
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by the magnetic diffusivity present in the simulations. 



We have quantified the scale-dependent anisotropy of the fluid motions in our convection 
simulations via measurements of the second-order structure functions. The flow is anisotropic 
at scales larger than 50 km and this appears to play a role in the large-scale magnetic energy 
spectrum generated by the dynamo in that it differs from the isotropic Kazantsev result 
and varies with the degree of anisotropy of the flow. The slope is significantly less steep 
than the k +3 ^ 2 law of the Kazantsev incompressible small-scale dynamo, but steepens with 
decreasing anisotropy of the flow as measured by the inertial-range velocity increments. This 
anisotropy decreases with Reynolds number as stronger horizontal fluctuations are generated 
by the turbulence (strong vertical fluctuations being induced by the convective driving). 

Our analysis suggests that solar surface convection is capable, via a turbulent small-scale 
dynamo, of generating and sustaining small-scale magnetic field. Of course, dynamo action in 
the Sun will also be present in lower layers where more kinetic energy is available and where 
stratification and rotation may also lead to an inverse cascade to large scales (a— effect). 
These two types of turbulent dynamos are likely intimately related. Even without a large- 
scale seed field local small-scale dynamo action should occur (e.g., during solar grand minima 
and for slowly rotating stars). In our simulation no large-scale background field is present. 
Turbulent induction of such a field (from a global dynamo) is an additional mechanism for 
the generation of small-scale field: it might obscure the presence of a small-scale dynamo. 
On the other hand, the large-scale magnetic field and the small-scale field produced from it 
will be subject to Ohmic decay: small-scale dynamo action might sustain magnetic field that 
would otherwise be lost. A future study should quantify at what background field strength 
the small-sc ale dynamo no longer dominates the production of small-scale magnetic field (as 
was done by lCattaneo et al.ll2003l for the Boussinesq case). 



The Pm <C 1 case will remain inaccessible for a long time to come for "realistic" small- 
scale dynamo simulations of the solar surface. Nonetheless, evidence of the existence of 



Pm <C 1 dynamos for idealistic simulations has now been given fjPonty et al.ll2005l ). When 
possible, radiative MHD magneto-convection calculations for Pm < 0.1 should be carried 
out to determine if the magnetic field generated has a different morphology for this case. 
Differences in the magnetic structuring has been demonstrated f or idealistic simulations o f 
small-scale dynamo action with Pm ^0.1 (jlskakov et al.l 120071 : ISchekochihin et al.l 120071 ). 
but these differences likely only affect the solar magnetic field at sub-kilometer scales. 

Small-scale dynamo action is no t suppressed for Pm -C 1 in idealistic simulations 
JPontv et alJl2005l : ; llskakov et all [20071 : ISchekochihin et al.ll2007h . We have shown that "re- 
alistic" simulations of dynamo action in strongly stratified, ra diative MHD with partia l 
ionization and little recirculation for Pm > 1 (first achieved by IVogler fc Schiisslerl 120071 ) 
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also do not suppress the small-scale dynamo nor do they supplant it with another mecha- 
nism. It is reasonable, then, to infer that small-scale dynamo action can occur for Pm <C 1 
combined with solar-like stratification, radiation, ionization, and recirculation. 
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A. Appendix 



A.l. Derivation of transfer functions 



Transfer analysis for the i ncompressible case is well-known (see iKraichnanl 119671 or for 
a more recent exposition, e.g., lAlexakis et al.ll2005l ). Here, we generalize the theory to the 
compressible and strongly stratified case. The latter requires a departure from the usual 
isotropic (and periodic) Fourier basis. Given a complete orthonormal basis, </>fc(x), and a 
function g(x), 

0(x) = £)0(AO&(x) (Al) 



with 



g(x.)4> k (x.)dx 3 



(A2) 



where Q is the analysis volume. Integrating the spatial density, g(x), of a global quantity, 
G, over the volume allows us to identify g(k) as a "spectral density" of G in k space, 

G= f g(x)dx 3 = [ J29(k)M*)dx 3 = J29(k) [ fc (x)dx 3 = V g(k) . (A3) 
Jn Ju k k Jn k 

For example, the magnetic spectral energy density satisfies 

E M = / e M dx 3 = , 



(A4) 



where bm = g^|B| 2 is the magnetic energy density (in Gaussian units) and E M is the total 
magnetic energy. From Parseval's theorem we see that 



f B(x) • B(x)dx 3 = V B(k) ■ B*(k) 



(A5) 
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where ■* specifies conjugate (depending on the basis). This allows us to identify that the 
magnetic spectral energy density is 

E M (k) = ±-B(k).B*(k). (A6) 

07T 

We can derive its temporal evolution (the magnetic component of the energy balance equa- 
tion) by projecting the induction equation, 

d t B + v VB = B • Vv - BV ■ v + r/V 2 B , (A7) 

onto the basis functions and then taking the dot product of this expression with -^B*(k). We 
add the conjugate of the result and divide the sum by two to derive that the time evolution 
of the spectral magnetic energy density, 

j t E M (k) = T KB (k)+T IB (k), (A8) 

is given by the magnetic energy transfer functions (representing the transfer of magnetic 
energy in k— space). In general, the transfer functions, Txy(k), measure the rate of energy 
transfer from energy reservoir X to the k— component of reservoir Y (also the work done 
on field Y by field X at wavenumber k). For Txy{k) > (< 0) energy is received (lost) by 
the k— component of Y. In particular, Tx B (k) = TxBT(k) + Tx B p{k) denotes the energy 
transfer rate from the kinetic energy reservoir to the k— component of the magnetic energy 
reservoir through the dynamo term (stretching against the magnetic tension force) J§ 

T KBT (k) = j-B(A;) ■ [B^W]*(fc) , (A9) 

and compression against magnetic pressure, 

T KBP {k) = -j-B(A;) • [BTv]>) - J-B(fc) ■ [v^VB]*(fc) . (A10) 

Tj B {k) represents the (negative of) energy loss rate to Joule heating, 

T IB (k) = ^B(k)-\y*B]\k). (All) 

47T 

The kinetic spectral energy density can be expressed 

E K {k) = ~ (v(Ar) • + **(*)) • (A12) 



6 Abbreviated expressions are used where the conjugate is assumed. The full expression is TKBr(k) 
± (B(fc) • [B^Vv]*(fc) + B*(fc) • [B^Vv](fc)). 
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Note that this is different from the velocity spectrum, 



E v (k) = -m-r(k) 



(A13) 



in the incompressible theory of Kolmogorov (see, e.g., lFrischlll995l ). The time evolution of 
the kinetic spectral energy density is given by 



jE K {k) = U{k) ■ d t \M\k) + \\M\k) ■ d t v(k) 



(A14) 



This is an alternative expression to the one derived by iMiura fc Kidal (119951 ) for the non- 
magnetic case. The kinetic energy transfer functions can be derived similarly to the expres- 
sion for the magnetic tr ansfer utilizing the conservation of momentum density equation, (see, 



e.g. 



VoglexetaU 120051), 



(A15) 



dt(pv) + V • (pvv) = -VP + pF + pV • r + -(J x B) , 

_ c 

and the conservation of mass density, p, 

d t p+V- (pv) = , 

which can be combined to form the momentum equation in terms of the velocity 

d t v + v • Vv = --VP + F + —V ■ t + —(J x B) . (A17) 
p p - pc 



(A16) 



Taking the appropriate dot products and integrating over the total volume, we derive that 
the kinetic part of the energy balance equation, 



d 
dt 



E K {k) = f(fc) + T KK {k) + T BK {k) + T IK {k) 



is given by the energy injection into the volume by body forces 

f(k) 



±v(k)-[ P F](k) + ±[pv]\k)-F(k) 



(A18) 



(A19) 



and the sum of the kinetic energy transfer functions. TxK{k) = TxKA{k) + TxKc{k) is the 
transfer of kinetic energy inside the kinetic energy reservoir ("cascade"), by advection 



T K KA(k) = 
and compressible (irrotational) motions 

T K Kc(k) = 



iv(*) • [v ^v)]\k) - l\pv](k) ■ [v^VvP) 



--v(A;).[pvV-v] (A:). 



(A20) 



(A21) 
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TiK:(k) = TjKc{k) +TjKD{k) represents energy transferred from the internal energy reservoir 
to the kinetic energy reservoir by compression, 



1 



1 



T IK c(k) = --v(fc) ■ [VP] (fc) 



-VP 



(*) 



(A22) 



and viscous dissipation, 



= --v(*) ■ /iV ■ r (k) - -\pv] (k) 



V-T 



(*)■ 



(A23) 



Tsxik) = TBKr(k) + TsKp(k) represents energy transferred from the magnetic energy reser- 
voir to the kinetic energy reservoir by the work of the Lorenz force both via magnetic tension, 



T BK T(k) = • pB^VB]*(fc) + ^[pvf(k) 



-B VB 

P 



(k) 



(A24) 



and via magnetic pressure, 



T B Kp{k) = 3-v(fc) 

07T 



■ 5 V|B|» 



W + (*) 



— VIBI 2 

2p 1 



(*). 



(A25) 



For an ideal gas with adiabatic index, 7, the internal energy density is given by e/ = 
P/(7 — 1) and the pressure is given by P = p9\T where T is the temperature and 9\ the gas 
constant. The internal spectral energy density can then be expressed 



91 



f(k)p*(k)+p(k)f*(k)) . 



2(7-1 

The time evolution of the internal spectral energy density is then 



(A26) 



d_ 
dt 



Ei{k) 



91 



7-1 



f(k)d t p*(k) 



91 



7-1 



P*(k)d t f(k) 



(A27) 



Using the conservation of mass, Equation (1A16I) . and the energy equation for temperature, 



91 



7 



- (d t T + v • VP) = V ■ (kVT) - PV • v + Q rad + p(r ■ V) 1 v + 



:f) 2± |J| 2 , (A28) 

47T (J 



we find the that time evolution of the internal spectral energy density, 

jE I {k) = U{k)+T I {k) 1 



(A29) 
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is given by the radiative heating (cooling), fi(k) = J n Qraddx 3 , and the sum of the inter- 
nal energy transfer function, Tj{k). Tj(k) = Tu{k) + TKi{k) + T B i(k). Tjj(fc) is the transfer 
of internal energy inside the internal energy reservoir: Tu(k) = — ^T(k)[V ■ (pv)] (A;) — 



T p*(k)[vVT](k) + p*(k) 



-V • (kVT) (k). T KI {k) = T KIC {k) +T KID (k) represents trans- 



fer of kinetic energy to internal energy by work against pressure gradients, Txic{k) 
-p*(k) 



iPV- V 

P 



(k) and viscous dissipation, T KID {k) = p*(k) ^p(r • V) T v (k). T BI {k) rep- 



resents transfer of magnetic energy to internal energy by Joule heating, T B i(k) = p*(k) ^{^) 2 ^\ 
where 



(*) 



J2T IB (k) +T BI (k) = -(^)- fv-(JxB)^ 3 . 
u An a J n 



(A30) 



As this is the integral over a divergence term, we recognize that the magnetic energy lost and 
the Joule heat gained differ only by a surface term. The sum of transfers between kinetic 
and magnetic energies is 



J2T KB (k)+T BK (k) = ^- ( V-((vxB)xB)^Efe 



(A31) 



where tyon is the inductive component of the Poynting flux into the domain. From this, we 
see that the time derivative of the total magnetic energy, 



Jt EM 



J2 T B(k) = Pan -J2 T BK(k) - ^T B/ (fc) 



(A32) 



is given by the inflow of electromagnetic energy (negative Poynting flux) at the boundary, 



Pan = -^- / V ■ (E x B) dx 3 
4?r Jn 



(E x B) ■ dA . 

47T 



(A33) 



the Lorentz force work against internal fluid motions, ^2 k T BK (k), and diffusive losses to 
Joule heating, ^2 k T B i(k). 

The right hand side of Equation (IA31I) is equivalent to a surface term via Gauss's 
divergence theorem, 

^ m = ^f ((vxB)xB)-dA. (A34) 
4vr J dn 

This surface term is zero for common choices of boundary conditions, hence, 



(A35) 
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This follows from the fact that Tkb and Tbk measure the transfer of energy between the 
magnetic and kinetic energy reservoirs. We can also show that 

y j T KB p{k)+T BKP {k) = ^- [ |B| 2 vdA. (A36) 

This is another surface term and 

J2 T KBp(k) = -J2 T BKp(k) (A37) 

k k 

as these two functions measure energy transfer via the interaction of compression and mag- 
netic pressure. In incompressible MHD, the second term of Tkbp in Equation flAlOj) . that 

T BB (k) = -— B{k) ■ [v • VB] (A;) , (A38) 

47T 

is identified as the transfer rate of magnetic energy to other scales within the magnetic energy 
reservoir (the magnetic energy "cascade"). As the first term of Tkbp is identically zero and 
^2k^BKp(k) = for incompressible flow, Ylk^Bsik) = (transfer is internal to magnetic 
energy reservoir). This is not true for compressible flow and the magnetic "cascade" cannot 
be separated from compressive kinetic energy transfers. From Equations (1A35j) and (1A37j) . 
we also find 

Tbkt (k) , (A39) 

k k 

as that these two functions measure energy transfer via the stretching of magnetic field lines 
against the magnetic tension force. A similar analysis shows, ^2 k TjB(k) = —J2k^Bi(k), 
J2k T iKD(k) = -J2k T KiD(k), and *E k T IKC (k) = -J2k T Kic(k). 



A. 2. Application to MURaM simulations 

For stratified convection, we use as the orthonormal basis functions a two-dimensional 
Fourier basis multiplied by vertical cardinal basis functions, 

to.,*,,* = ^e-^ x+k y^6(z - z') . (A40) 

To study the generation (or loss of) magnetic energy at a given scale, the resulting transfer 
functions are projected onto a one-dimensional wavenumber, k\ = k^ + ky, and summed over 
the vertical direction, 



T K B(k±) = ^2 ^2 T KB (k x ,k y ,z) 



(A41) 
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To quantify the source of magnetic energy generation at the scale 27r/k±, we need only 
measure T^bt and T^bp- The resistive transfer, Tib, mostly serves as a sink for magnetic 
energy at small (< 5Ax) scales. To complete the physical picture, we also measure T BKT and 
Tbkp-i the rate at which kinetic energy is lost to (gained from) the magnetic field. Transfers 
between internal and kinetic energies may be ignored as we are only interested in generation 
and loss of magnetic energy in a dynamo analysis. As shown in the previous section, Equation 
(IA3ip . the sum of these four transfers is equal to the inflow of electromagnetic energy from 
the inductive component of the Poynting flux, tygn, into the domain, 

Trbt + T KBP + T B kt + Tbkp + T MM + A chain = ^Pqq . (A42) 

Zbot< z < z top fc± 

Here Tmm(^) = ■ [vV • B] (A;) denotes the influence of the numerical artifacts of 

magnetic monopoles and 



*chain 



v • Ve M - • (v ■ VB) dx 3 (A43) 

47T 



n 

2 



arises because of numerical inaccuracies in the chain rule, 0.5dB /dx = BdB/dx. Multiply- 
ing MURaM's 5-pt stencil for the derivative of B by B yields, 

dB _ Aa^ (8) (A44) 

dx 30 v ; 



B ^± - ^-{Q^-B^ + 12^B^ + 5B^B^) . (A45) 



while the numerical derivative of B 2 /2 yields 

-^-^-(6^ 4 > + 12^ 
dx 36 dx dx 2 

As strong gradients of the magnetic field exist in the downflows, o(B) -C o(B^) and sig- 
nificant departures, ~ 8% of the total transfer rate, from the chain rule occur. We track 
this through direct calculation of A c h a in- Lower order errors occur for the two grid points 
closest to the top and bottom boundaries. For this reason we choose our analysis domain 
such that the boundary, dQ, is halfway between the second and third grid points from the 
bottom (and top), tydn is interpolated using Newton divided differences and the error in 
calculating Equation (IA42j) is computed for each snapshot. As Equation (1A43[) accounts for 



the chain-rule error only for the largest and most easily calculated contribution, some small 
residual error is found. This error, always less than 3%, is assumed to be spread equally over 
k±— space and is reported for all plots of transfer functions. 
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